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A Nonequilibrium Model for a Moderate Pressure Hydrogen Microwave Discharge Plasma 

Carl D. Scott* 

Laboratoire d'Ingenierie des Matdriaux et des Hautes-Pressions 
Universite Paris Nord 
Av. Jean-Baptiste Clement 
93430 Villetaneuse 
France 


ABSTRACT 

This document describes a simple nonequilibrium energy exchange and chemical reaction model to 
be used in a computational fluid dynamics calculation for a hydrogen plasma excited by 
microwaves. The model takes into account the exchange between the electrons and excited states of 
molecular and atomic hydrogen. Specifically, electron-translation, electron-vibration, translation- 
vibration, ionization and dissociation are included. The model assumes three temperatures, 
translational/rotational, vibrational, and electron, each describing a Boltzmann distribution for its 
respective energy mode. The energy from the microwave source is coupled to the energy equation 
via a source term that depends on an effective electric field which must be calculated outside the 
present model. This electric field must be found by coupling the results of the fluid dynamics and 
kinetics solution with a solution to Maxwell's equations that includes the effects of the plasma 
permittivity. The solution to Maxwell's equations is not within the scope of this present paper. 


NOMENCLATURE 

Cj mass fraction 

D ls binary diffusion coefficient 

£ electric field vector 

e electron charge 

e total internal energy per unit mass 

e*(T) equilibrium vibrational energy per unit mass evaluated 
at gas temperature 

Ee energy per unit volume of electrons 

e e energy per unit mass of electrons 

Ejf ionization energy for the f* species 

Eq threshold energy for a reaction 

E v vibrational energy per unit volume 

ev vibrational energy per unit mass 

f s body force on species s 

h^ enthalpy of formation of species s 

j vibrational quantum number 

kb recombination rate coefficient 

k 0 j rate coefficient for excitation of jth energy level by electrons 

kim rate coefficient for ionization of molecules by electrons 

kia rate coefficient for ionization of atoms by electrons 

m mass of particle 


* This work was done while Visiting Professor. Permanent address: EG3, NASA Johnson Space Center, Houston, 
Texas 77058. 
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M s 

nef 

n s 




Qd 

4 

4 

Qe 

q e 

Qe-v 

Qrad 

Qr-e 

Qt-v 

Qt-v 

a 

q v 

T 

T e 

T v 

u 

V e 

V s 

W 

W e 

Ws 

xi 

X s 

Z s 

Greek 

Sij 

y 

x 

Xe 


mass flux vector of species s 
number density of species f 
number density of species s 
pressure 
electron pressure 

energy exchange rate per unit volume from electrons to excited 
singlet states 

energy exchange rate per unit volume from electrons to 

vibrational ground electronic states via singlet state 
excitation 

electron energy exchange rate per unit volume due to 
dissociation 

electron energy loss rate per unit volume for excitation of 
triplet states 

rate of exchange of energy to electrons per unit volume 

rate of energy exchange per unit volume to vibration 

rate of energy gain per unit volume from the electromagnetic field 
electron energy flux vector 

electron-vibration energy exchange rate per unit volume 
radiative energy exchange with gas (usually a loss) 
electron-translation energy exchange rate per unit volume 
translation-vibration energy exchange rate per unit volume 
for molecule-molecule collisions 
translation- vibration energy exchange rate per unit volume 

for atom-molecule collisions 
vibrational energy flux vector 
translational/rotational temperature 
electron temperature 
vibrational temperature 
gas velocity 

electron diffusion velocity 
diffusion velocity of species s 
molecular weight of mixture 
molecular weight of electron 
molecular weight of species s 
jth component of position vector 
mole fraction of species s 
degree of ionization of species s 


Kronecker delta 

wall recombination probability 

thermal conductivity 

thermal conductivity for electrons 
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X' e thermal conductivity for electrons considering electrons only 
Xy thermal conductivity for vibrational energy 

X'v thermal conductivity for vibrational energy considering only vibrating molecules 
(i viscosity 

Vc(u) electron-heavy particle collision frequency at energy u 
v~ heavy particle-electron momentum transfer collision frequency 

v^ heavy particle-electron momentum transfer collision frequency 

weighted with electron energy 
6 energy loss rate per electron 

p mass density of gas 

Ps mass density of species s 

0 cross section 

<% electron-heavy particle momentum exchange cross section 

Of real part of complex electrical conductivity 

x shear stress tensor 

Te-v electron-vibration relaxation time 

tv translation-vibrational relaxation time 

x v translation-vibrational relaxation time H 2 -H 2 

Hj-Hj 

Tyjs vibrational relaxation time for pair 1 and s 

to frequency of microwave 

cos rate of production per unit volume of species s 

Subscripts and Superscripts 

e electron 

1 ionization 

r heavy particles 

s species s 

T translational 

v vibrational 

INTRODUCTION 

Hydrogen plasmas with a few percent of methane and a trace of oxygen have been used to produce 
diamond films in the laboratory. The rates of deposition and etching of graphite and diamond on 
substrates are significantly affected by the concentrations of the various species and the energy 
states in which they exist Likewise, the surface temperature and flow pattern in the reactor affect- 
transport of reactive species to the surface. Therefore, a program to compute the chemistry, flow 
field, and energy states of a hydrogen plasma has been started. Eventually, the flow simulation 
must include the carbon containing species and oxygen, but for the present we have started with a 
pure hydrogen microwave generated plasma to be discussed in the present paper. 
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The subject of modeling the energy exchange between the various energy states in hydrogen 
plasmas has been developed over die last two or three decades with various degrees of complexity 
depending on the plasma conditions. Cross sections and reactions rates for electron-hydrogen 
collisions have been measured and calculated by a number of workers. These interactions involve 
dissociation, excitation, electron- vibration energy exchange, ionization and translational-vibrational 
energy exchange. The various paths of energy loss from electrons to the internal degrees of 
freedom of hydrogen and then to the translational/rotational degrees of freedom affect the ultimate 
conditions of the plasma and thus the possible chemical reactions both in the gas phase and on 
surfaces. 

This paper defines equations and reactions for modeling experiments 1 atthe Laboratoire dlng£ni&ie 
des Mat6riaux et des Hautes-Pressions at the Univerist6 Paris Nord, Villetaneuse, France. In these 
experiments a pressure level of about 4 kPa assures a reasonable rate of deposition. The plasma is 
established by a microwave system operating at about 500 watts. A schematic of the apparatus is 
given in Fig. 1. The plasma zone that is approximately hemispherical in shape exists next to the 
substrate surface. This surface is essentially immersed in a region of high electric field produced by 
a focused microwave antenna. The size and conditions of the plasma depend on the power coupled 
to this cavity as well as the electron and energy loss mechanisms in the plasma. These loss 
mechanisms include diffusion of electrons and excited species to the surface and out of the high 
microwave field zone, electron energy loss by collision with H 2 to form vibrationally or 
electronically excited molecules, dissociation of the molecules, ionization of atoms and molecules, 
and electron-translation/rotation exchange. Forced and free convection also play a role in loss of 
electrons and other excited species from the plasma zone by removing electrons and other excited 
species from the active plasma zone. 

Modeling of hydrogen plasmas and the study of energy exchange mechanisms have been done by 
several workers in other applications. Among them are the works of Gorse, et al. 2 , Gorse et al., 3 
who applied their models to multicusp magnetic discharge plasma, and Loureiro and Ferreira, 4 who 
applied their model to the positive column of a dc discharge. These applications are at sufficiently 
low temperature and pressure that it was necessary to account for the strong nonequilibrium 
associated with the electronic energy distribution function. Their distributions were expected to 
deviate from a Boltzmann distribution. In the present case the pressure and temperature are 
sufficiently high that the electron energy distribution function may be nearly Boltzmann, and the 
vibrational distribution function of H 2 can expected to be Boltzmann also. However, it is expected 
that the temperatures describing these Boltzmann distributions for the various degrees of freedom 
are not the same, nor are they equal to the translational/rotational temperature. Validation of these 
expectations will await spectroscopic and other diagnostic measurements. 

The model will be introduced into a finite difference solution of the Navier-Stokes equations set up 
for the geometry of the plasma reactor. In lieu of a solution of measurements of the electron 
temperature and density. It will also be necessary to couple the microwave power into the model via 
results from a solution of Maxwell's equations. This procedure must involve an iteration with the 
solution of the Navier-Stokes equations because there is a strong coupling between the loss 
mechanisms computed by the N-S solution and the energy coupled to the plasma from the 
microwave source. The electromagnetic field of the source and in the plasma depend strongly on 
plasma conditions. 
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HYDROGEN PLASMA ENERGY EXCHANGE MODEL 
The Flow Equations 

It is assumed that the microwave discharge plasma and the associated flow in the reactor can be 
described by the continuum Navier-Stokes equations in a form patterned after those presented for 
shock heated air by Lee 5 . Since there are several modes of energy that are not strongly coupled we 
expect that it is necessary to account for three temperatures and thus use three energy equations. 

The vibrational distribution of the molecules is assumed to be describable by a vibrational 
temperature, the electrons and electronic excited states are assumed to be described by an electron- 
electronic temperature, (called electron temperature in subsequent paragraphs), and the heavy 
particle translation and rotation distribution can be described by a translational! rotational 
temperature. ’Hie chemical species are not assumed to be in equilibrium, but must be described by 
chemical kinetic relations. In addition to the momentum equations and continuity equation, whose 
forms are independent of the kinetic model, we have the set of species equations, one for each 
species . The three energy equations are electron, vibration, and total. These equations whose form 
is dependent on the kinetics, are given here. The time dependent form is given even though the 
problem is actually steady state because the solution technique may be a time relaxation scheme in 
which the parabolic coordinate is time. Otherwise, the steady state equations are elliptic, which 
makes the solution technique for the problem less flexible. 

Species equation: 


+ V-(p s u + PsV s ) = co s 


Electron energy equation: 

Tt 6 + v (E e u) = -V q c + Q* x - V ((Ec + p e )V c ) 

- V-(peti) - enetic f +Q e 

Vibrational energy equation: 


( 1 ) 


( 2 ) 



+ V(E v u) = -Vq v + Q* x 


n 

^■Xp* e vV, 


i=i 


(3) 


Total energy equation: 


■Jf + V*(peu) = -V q -V-(up) + V (u x) + Xp s u s f s + Q e - Qmd 

s=l 


(4) 


where, p s is the density of species s, u is the mass averaged convection velocity. p s V s is the 
diffusion mass flux of species s given by 


PsV,= P^2> u [VC, + C|Vln(W)] 


(5) 
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where W s and W are the molecular weights of species s and the mixture, respectively, and is 
the multicomponent diffusion coefficient for the l,s pair of species, and Q is the species mass 

fraction pi/p. 

The rate of production of species s is coj. Ee and Ey and q e and q v are the specific energies and 
energy fluxes of the electrons and vibrating molecules respectively. The total internal energy per 
unit mass is e. The energy per unit mass of electrons and of vibrational energy are denoted ee and 
ey respectively. The rates of exchange of energy per unit volume between the various modes and 
vibration is given by 


Qgj = Qt-V + Qt-v^ + C^; v + Qe-v (6) 

and the rate of exchange to the electrons is given by 

= Qr< -Qe-v - 0,! e v - E i f < 7) 

which exhibit the rate of exchange of energy per unit volume between the various modes designated 
by die subscripts. The energy loss by the electrons in exciting the singlet states is <^ v and the 

lv 

energy gain by the vibrational states in this process is C^y and their difference accounts for the 
radiation loss from the singlet states of H 2 and contributes to Qrad- The energy lost by electron in 
exciting the triplet electronic levels which lead to dissociation is <^j. The last term in (7) is the rate 

of energy loss per unit volume by electrons ionizing the f* species, in the present case forming Hj 
and H + . 

The shear stress tensor is denoted X 



The external forces are denoted f s . 

The viscosity and bulk viscosity are respectively |i and )ib> The heat flux vectors can be 
expressed by 


q — -XVT -X*VT e -XyVTy + Sps^s (9) 

q c = -X’cVTe (10) 

qv = -XyVTy 


( 11 ) 
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The thermal conductivities of the gas, the electrons, and the vibrational energy are, respectively A, 
A© and A v , while the modified vibrational and electron thermal conductivities A' v and A' e are 
defined by Lee 5 . 

Electrostatic forces are the only ones that may be present due to induced electric fields. This force 
can be expressed by 




p s f s ~ n s eZs£ 


( 12 ) 


where £ is the time averaged induced electric field, which it is likely to be zero in the continuum 
regime for usual conditions of microwave plasmas, except in regions of strong gradients. The 
electric field can be computed by the Langmuir and Tonks formula 


£ 



( 13 ) 


The microwave source pumps energy into the electrons through an electric-magnetic wave that 
penetrates the plasma. This e-m heating depends on the time averaged electric field which depends 
on the plasma conditions as well as the electromagnetic and geometrical characteristics of the source. 
Thus as mentioned in the introduction the fluids problem is coupled to the electromagnetic problem. 
In the present treatment of the plasma it will be necessary to solve (he problem iteratively by 
calculating Q e (x) by assuming a distribution of the mean electric field in the plasma zone, then 
updating its value by solving Maxwell's equations with the resultant plasma conditions from the 
solution to the fluids equations. An alternative to solving Maxwell's equations is to use measured 
electron temperature and density distributions when computing Qe(x). However, measuring these 
quantities is difficult, particularly without resorting to Langmuir probes which could perturb the 
plasma. 


Modeling Energy Exchange Between Modes 

Energy exchange between the electrons, vibration, translation/rotation, and dissociation is accounted 
for by the terms Qr-V» Qe-V> Qt-c. Qe-V» and Qe-d. which will be defined in this section. 

Translation- Vibration Rel axation (Molecules') . 

One could calculate the T-v exchange for molecules by summing over all rate probabilities as shown 
in Cacciatore, et al. 6 , and basing the rates on the data of Billing and Fisher 7 and the method of 
Audibert, et al.,® but we shall adopt some of the modeling methods of Lee, 5 who applied Landau- 
Teller type relations for vibration-translation relaxation to mixtures of nitrogen and oxygen. The 
energy exchange rate is given by 


Qt-v 


(P HJ + P HJ +) 


<(T> - e. 


Tv 


f <(T) - e v 

* rr\ 

^e v (T) 0 - Cyo 


s-1 


( 14 ) 


where e v (D is the equilibrium specific vibrational energy at the translational temperature, ey is the 
vibrational energy of the flow, and x v is the relaxation time for T-v relaxation. The subscript o 
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indicates the initial value (t = 0) of the energy and s £ 1 is an exponent that takes the value 
determined by a fit of the relaxation with exact solutions of the relaxation rate 9 . Lee 5 has given the 
general form for relaxation of species 1 in a mixture. 


T V 1 - 



(15) 


The t v , values can be determined from cither the Millican and White 10 or Kiefer and Lutz 11 
Is 

correlations. Audibert, et al.6 have shown that the Kiefer and Lutz 11 correlations agree with their 
higher temperature results over a fairly broad temperature range of interest here. The Kiefer and 
Lutz 11 correlation is 


t v p = 3.9xlO* 10 exp[Y^j atm-s 


(16) 


For now we use only the one expression since there are no other molecules in the system. The rate 
of energy exchange between translational and vibrational modes is 


3e c v (T) - Cv 

Qr-v - Ph 2 !h “ p h 2 

V «2-H2 

where e*(T) is the equilibrium vibrational energy per unit mass evaluated at the translational 
temperature and e v is the vibrational energy per unit mass. 

Translation-Vibration Relaxation (Atoms) 

The relaxation scheme for atoms exchanging translational energy with molecules involves 
significant quenching of hydrogen molecule vibrational states. The rates are given in the appendix 

of Gorse, et al. 3 The relaxation expression is of the same form as equation (14) and the two terms 
are additive. There are two principle mechanisms that transfer energy from vibration to translation 
when H atoms collide with H 2 molecules. The first is an echange of energy without any chemical 
rearrangement of the molecule. The other involves the exchange of an atom by first forming a 
complex of three atoms, which immediately split apart, transferring vibrational energy in the process. 
Rates for these two reaction paths are given in Ref. 3. There are rates for many combinations of 
initial and final vibrational energy, but since the ground state dominates at low temperatures we only 
consider transistions involving the ground state and the first excited vibrational state of H 2 . We 
have included the transition 1,0 for the two paths in Table 2. The energy transfer expression for 
this exchange is given by 

Qr-Vfi — (k”r-v,a + hT-v^r) n H 2 n H (1®) 

where kr-va and kj-v^r are the direct and exchange energy transfer rates respectively and Ey is the 
energy of the v=l state. If conditions are such that additional vibrational states are important, these 


rates can be easily added, weighting them with the relative populations determined by the vibrational 
temperature and summing. 

Electron-Translation Exchange 

The form of the electron-translation energy transfer rate Qr-e is taken from Lee-* with relaxation 
coefficients calculated from the momentum exchange cross section data of Buckman and Phelps. 12 
The momentum exchange collision frequency is determined from integrals of the cross sections over 
the electron velocity distribution function 


** 


1 JtferC'efe ^ c 'e 


or the energy averaged collision frequency 


• _ “em- 
er 3neia { 


/' 


tferC-fe dc' c 


(19) 


( 20 ) 


These collision frequencies, written in terms of a Boltzmann energy distribution function, can be 
simplified to 


v‘*n/s] 

a 6.69x1 0 7 


and 


=TrppJo«{om 2 ]E[eV] e xp(.BTJdE 


v^ri/si 


^3J - K[cm2jE[eVl2cxp(.™ dE 


( 21 ) 


( 22 ) 


Hiese collision frequencies are plotted against electron temperature in Fig. 2. Polynomial curve fits 
of the results are given in the figure and in Table 2. 

The elastic collision transfer term Qr-e is given by Sutton and Sherman 1 ^ by the relation 




Or* “*^( 1 * J 0c,c e f e dc ’ 

- mcn r (u c - u r ) u g Ja er c' c f e dc’ c 
Lee 7 expressed this in terms of the collision frequencies as 


(23) 


Qr-c = 2mene 


?kT c 


) 


(24) 


where a term involving the relative velocity (diffusion) between electrons and heavy particles has 
been neglected because of the assumption of zero current in the plasma. For a simple gas Kir* 
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hydrogen only two species, H 2 and H are considered because the ionization level is very small. 
Since the collision cross sections for momentum exchange were obtained for a mixture of H 2 and 
H, an average cross section was obtained in the determination of the cross sections, thus we will 
only consider equation (24) as having one term with nv equal to the mean mass. 

Electron- Vi bration Exchange. 14 

As indicated in Table 1 there are two paths to vibrational exchange with electrons. The first is direct 
vibration excitation and deexcitation in a given electronic state — in particular, the ground state. The 
second is vibrational excitation via electronic excitation of excited electronic singlet states followed 
subsequently by a radiative transition to a vibrationally excited ground state. Cross sections for 
these reactions are given in Ref. 12. The rate coefficients for these reactions are calculated by 
integrating the cross sections over the electron energy distribution function. 


koj(T e )[cm3/s] = 6.70x107 t; 3/2 Ja 0 j(E)Eexp(-E/T c )dE (25) 

Eo 

where the energy E and the temperature are in electron- Volts. Eq is the threshold energy of the 
exchange, and a 0 j is in cm 2 . 


The relaxation time Xe- V is determined following the procedure developed in the appendix of Ref. 5. 


Te-v — 


1 

ne(l - exp(-e v /T c )) 2 |Jk 0 j j 2 dj 


(26) 


It has been assumed that the vibrational levels are closely spaced so the integral in (26) well 
represents a sum over levels. This may not be the case for H 2 , but we will use it until the error is 
determined to be too large. The form of kqj for the first three vibrational levels (j=l-3) are very 
similar in shape. Therefore, an approximation was made that the integral in (26) could be done 
semi-analytically by fitting the maxima of the koj to a simple function of the form 

koj max = 3.39x10-8 e- 2 - 401 i (27) 

which can be integrated analytically. The peak values are shown in Fig. 3. A polynomial curve fit 
of the result is shown in the figure. The relative rate coefficient function was integrated numerically 

to obtain a constant factor. The resulting relaxation time Xe-V as function of electron energy is 
given in Fig. 4. The loss rate of energy is obtained from the relaxation time by the relation 


dev 


Qc ' v = Ph 2 ’^ = Ph : 


c y (T c ) - e v 
W 


r * 


v 


c v (T) - ey 

Cy(T)o - Cvo 




s-1 


J 


(28) 


where e (T c ) is the equilibrium vibrational energy evaluated at the electron temperature (the 
driving energy) and e v is the current vibrational energy. 

The excitation of vibrational energy via the second process (process no. 3 in the Table 1) that 
involves the singlet states leads to a loss of energy by the electrons since the plasma is optically thin. 


It results in a gain by the vibrational states of the ground electronic state. A negligible amount of 
this radiation is absorbed by the plasma and thus the process is hot reversed. 

The loss rate by the electrons is then given by 

= nen^Emkim (29) 

m 

The gain by the vibrational energy is similar, being 

v = n e n Z(Em'Emv)k i m (30) 

m 

where k 3 nj are the rate coefficients for the excitation of the ground electronic state of H 2 to the m* 
singlet excited state and E m is the energy per molecule required to excite the molecule and E mv is 
the mean vibrational energy of the final state to which the molecule radiates. The radiative energy 
loss due to this process is the difference between these terms 

QrxU'Q^V'Qe-V < 31 > 

Electronic state excitati on bv Electrons 

Elections may lose energy via exciation of electronic states if atoms or molecules. This has already 
be addressed for excitation of singlet and triplet states of H 2 which leads to vibrational excitation 
and dissociation respectively. Atomic hydrogen may also be excited by electron collisions. Rates 

for these reactions can be found from data in Janev, et al. 15 Excitation of the n=3 state of H via the 
reaction 


e + H -> e + H(n=3) (32) 

and 

e + H 2 e + H(n=3) + H (33) 

are given in Table 2. This leads to radiation from the n=3 to the n=2 state which should be included 
in the term Qrad- Due to resonance of the n=2 state, we will neglect radiation cascading from it, so 
that the radiative energy loss rate is given by 

Qrad,n=3 = n e(kn=3,a n H + kn=3,m n H2^3,2 


Energy Absorbed from Microwave Electromagnetic Field 

The energy absorbed by the plasma from the electromagnetic field can be obtained by the technique 
of Ferreira, Moisan, and Zakrzewski 16 . In a global sense the power absorbed is exactly equal to the 
energy lost by electrons to die other particles and to the walls of the container. Losses are due to 

electron collisions with the other species and the wall. In their book 15 the energy loss rate per 
election in the plasma is given by 


®L * ~ £ <Vc(E)E> + ^^<v n (E)>e^ >n + <Vi(E)>eIk. 


(34) 
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where the brackets "< >" denote averages over the electron energy distribution function. The 
threshold energies for excitation and ionization are I^ n and 1^. , respectively. The plasma will be 

in steady state when the power absorbed by the plasma n e 0A equals the power lost by the electrons 
ne0L. The local absorbed power is given by 

QeOO = He® A = ®r(£'Orms (35) 

where is the real part of the complex electric conductivity 


o r = 


2 e? f v c (E) 

3 me J v c 2 (E) + CO 2 

0 


"“If dE 


(36) 


F c is the local electron energy distribution function and v c (E) = na er (E)V2E/m e is the collision 
frequency for electrons of energy u. The relations given in equations (35) and (36) can be used to 
compute the source term Q e (x) in the energy equations given the mas value of the electric field. 
The electric field must be found from a solution of Maxwell's equations. If the electron density and 
temperature distributions are known then one can compute the local value of the energy loss rate 

6l> hence Q e (x) without having to solve Maxwell's equations. In that case one must assume that 

0L = 0a locally. However, this leaves the problem of boundaries where there may be loss of 
electrons by convection or diffusion. 


CHEMICAL REACTION MODEL 

The gas phase chemical reactions that are assumed in this model fall in two basic categories. First 
are those that involve reactions with electrons, e.g., dissociation and ionization, while the other 
category includes the recombination of ions and atoms. 

Electron and Ion Chemistry 

The reactions for electron and ion chemistry are summarized here as 


e + H 2 (v) e + H 2 * -» e + 2H 

kd 

(37) 

e + H 2 (v) -> e + H 2 * -4 e + H + H* 

kd* 

(38) 

e + H 2 (v) -» e + + e 

ka 

(39) 

c + H-»e + H + + e 

kii 

(40) 

2e + H+-*e + H 

krl 

(41) 

2e + H 2 — ♦ e + H 2 

kr2 

(42) 
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e + H + hv + H 

krrl 

(43) 

— 

2e + H^v) e + H + H 

kdrv 

(44) 


e + h£ -» hv + H 2 

krr2 

(45) 


The reaction rates for dissociation (37) and ionization (39) are calculated from cross section data 
from Ref. 12 by integration over the electron energy distribution function that is assumed to be 
Boltzmannian at an electron temperature T e . Excitation of the triplet states of H 2 , particularly the 

b 3 £j, c3n„, a3 -£g* ®nd d^n states, lead to dissociation. These reactions are denoted by reaction 

(37). Dissociation rate coefficients computed from the cross sections of Ref. 12 using equation (25) 
are given in Fig. 5, along with an Arrhenius curve fit. 

A plot and curve fit of the ionization rate coefficient for H 2 are given in Fig. 6, where the cross 
sections were obtained from Ref. 12. Although the cross sections were obtained and presented as 
functions of l£l/n where a Boltzmann distribution did not obtain, we have assumed that the 
Boltzmann distribution of electrons is the case here and we have used the cross sections of Ref. 12 
directly. 

For sufficiently energetic electrons and with H 2 molecule in excited states a dissociation as indicated 
by (38) can result in the production of an excited atom. Cross sections and rate coefficients for 
tiiese reactions have been calculated by Celiberto, ct al. 17 For the present, these reactions shall not 
be incorporated into the model. However, it may be necessary to account for these reactions for 
diagnostic purposes because the other path for excitation of H atoms, direct excitation of H atoms 
by electron impact, may not be the only significant path of excitation of the Ha and Hr lines, from 
which we may estimate the electron temperature. 

The reverse (recombination and deionization) reactions are denoted by reactions (41)-(44). The rate 
the rate of deionization is evaluated at the electron temperature T e . Recombination of hydrogen 
atomic ions can follow two paths, one is three body recombination as indicated by (42), whereas 
the other is radiative recombination as indicated in (43). The rate of three body recombination 
depends on the energy level in which the H-atom is found, thus on the kinetics of the problem and 
the electron density and temperature. For sufficiently high electron temperature and low electron 
density the recombination coefficient is fairly insensitive to the electron density because it is 

dominated by radiative recombination. Radiation trapping of the Lyman-a line can affect the rate of 
recombination since the number density of the ground and excited states of the atoms may be 
affected. Drawin and Emard 18 have solved the master equations, taking into account both three 
body recombination and radiative recombination with radiative absorption. Their results show a 
dependence on electron density in lower temperature and higher density regimes. Their results are 
shown in Fig. 7 along with the purely radiative recombination rate from Massey 19 and purely three 
body recombination rate from a curve fit of Johnson and Hinnov 7 ® measurements 71 for relatively 
high electron density. The path of molecular ion radiative recombination (45) is not included in the 
model because of lack of rate information. This path may be improbable compared with dissociative 
recombination. 


The molecular and atomic hydrogen ions may undergo an exchange reaction 
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H+ + H -> H 2 + H + (46) 

where the reaction rate for production of H + is given by Kaipas, et al. 22 in Table 3. 

The rate of recombination of with electrons to form a neutral molecule is not known (i. e., has 

not been found in the literature) which probably means that the encounter of electrons with 

results in dissociative recombination, reaction (44) instead of direct recombination. The rate for this 
reaction is from Nakashima, et al. 23 

Two reactions that should be included that involve the ions H* and H 2 + are 


e + H 2 -> H- + H(ls) 

kda 

(47) 

H 2 + + H 2 -»H3 + + H 

ke23 

(48) 


The first, (47) is dissociative attachment and it is significant in hydorgen plasmas. The other is the 
productionof the H3 + ion (48). The rates for these reactions derived from data in Janev, et al. 15 are 
given in Table 3. 

Neutral Gas Chemistry 

The only reactions in a low gas temperature pure hydrogen plasma that belong to the category of 
neutral gas chemistry are the recombination and dissociation reactions given by the equations 


2H + M — » H 2 + M 

kb 

(49) 

H 2 + H 2 -*H + H + H 2 

kd 

(50) 

h 2 + h-»h+h+h 

kd 

(51) 


The rates for this three body recombination are found in Cohen and Westburg's 24 data sheets, and 
the rate of the recombination reaction is evaluated at the gas temperature T. The rates of the 
dissociation reaction involve vibrational excitation and therefore depend on the vibrational 
temperature in a way not yet determined. 

The model proposed for computation of the hydrogen/electron chemistry and energy exchange in the 
mic rowave plasma is s ummar ized in Table 1. Expressions for collision frequencies and relaxation 
times to be used in energy exchange calculations are given in Table 2. Rate coefficients for the 
chemical reactions are given in Table 3. 

Gas/Surface Reactions 

The reactions that occur on surfaces in a hydrogen plasma depend to some extent on the surface 
materials It has usually been accepted that almost 100% of incident ions recombine catalytically on 
surfaces. We will assume that is true in the present model. The neutral hydrogen atoms may react 
catalytically on man y surfaces to form the hydrogen molecule. There is evidence that the molecules 
that are formed on me tals are in excited states. These states are probably not in thermal equilibrium 
with the gas or the surface and it may not be possible to ascribe a temperature to these states. It 
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remains an area of research to determine what effects these excited nonequilibrium species have on 
the flow and its chemistry. 

Another category of reactions that may occur on materials is chemical reaction with the surface 
atoms and molecules. These reactions may etch away material such as hydrogen or oxygen etching 
of graphite by forming chemical bonds and volatile products. They also may be involved in the 
formation of phases such as diamond when carbon is in the mixture. For the present pure hydrogen 
case we will assume that no chemical reactions involving hydrogen take place. 

The reactions assumed in the model are the following catalytic reactions involving ions and atoms. 
Excited state production and quenching shall be indicated in die reaction equations, but the reactions 
rates for them shall remain unknown for the present. 


H+ + S — >H + S (52) 

or possibly H+ + H-S H 2 (*,v) + S (53) 

H+ + S -> H 2 + S (54) 

H + H-S -> H 2 (*,v) + S (55) 

H 2 (v) + S -4 H 2 (w) + S (56) 

H 2 (v) + S -> H 2 (0) (57) 


The surface materials in the present experiments to be modeled are silicon, graphite, diamond, and 
adsorbed hydrogen. 

Wall Boundary Conditions for the Species Equations 

Hie wall boundary conditions for a low density partially ionized air were developed by Gupta, 
Scott, and Moss.^ The boundary conditions for the ions and atoms in an eleven species air model 
were subsequendy simplified by Scott^ for a nonslip model. These models are phenomenological 
in that they do not address chemical mechanisms of the reactions nor the physical processes taking 
place. One must know from some source the overall reaction probability as a function of wall 

temperature or other parameters. The probability of recombination y, is simp ly defined as the 
fraction of the incident atoms that recombine on the surface. This means that the reaction rate can be 
written as 


c w*Y 


kT 

2rcm 


and the net normal flux of atoms due to recombination is 

Ma = pCaY- 


V 2jcm 


(58) 


(59) 


for a first order reaction. The net flux of species i is written in terms of the diffusion of species as 
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Mj = pVj = p (60) 

where Dy is the multicomponent binary diffusion coefficient, Q is the species mass fraction, and p 
is the total mass density. By equating (59) and (60) we get a mixed boundary condition for species 
i and the other species j. The boundary condition is mixed in the sense that it contains both Q 
and derivatives of Q. If more than one surface reaction can produce or consume species i, then 
equation (59) will have additional terms that include these reactions. For example, the 
consumption/production of H atoms on a surface may involve the recombination of the H atoms and 
the recombination of the H + ion, viz. 


M H = -C H 7h 


V 


kT 


SL 


2jimH 


Sr 


v 


kT 


J8L 


27raiH 


(61) 


where we have set the ion recombination coefficient equal to unity. The first term on the right 
corresponds to consumption of hydrogen atoms by recombination, reaction (55), and the second 
corresponds to production of hydrogen atoms by hydrogen atomic ion recombination, reaction (52). 
If atomic hydrogen ions recombine via reaction (53) then equation (61) would have only the first 
term and the corresponding to equation for H 2 would include the atomic ion recombination term in 
addition to terms for atom recombination and molecular ion recombination, viz. 

^ - ch * V£T + V5T + V- <“> 

Note that (61) and (62) are not consistent with each other since they assume different paths for 
recombination of H+. They are given her as examples of possible species equation boundary 
conditions. 


CONCLUSIONS 

A model of a hydrogen plasma at moderate pressures has been proposed for use in calculating the 
properties of a microwave discharge used for developing diamond film deposition techniques. The 
model assumes thermal nonequilibrium between die rotational/translational, vibrational, and electron 
degrees of freedom, but assumes that each degree of freedom can be described by a Boltzmann 
distribution at a given temperature. Chemical nonequilibrium, including ions, is incorporated in the 
model. Simple catalytic gas/surface interactions are included. The model is fashioned after that of 
Lee for air at hypersonic flight conditions and includes specific hydrogen reaction and excitation 
rates. The model can be included in a flow simulation calculation such as a finite difference solution 
of die Navier-Stokes equations. Validation of the model must be based on its ability to predict 
measurements in the hydrogen microwave plasma. The source term for input of energy from the 
microwave source must be obtained either from a solution to Maxwell's equations for die same 
geometrical set up as the fluid dynamics equations, or must be obtained from measurements of die 
electron density and temperature. As confidence is gained for the model and the computational 
method for pure hydrogen, additional species will be added to permit one to account for deposition 
of carbon in the forms of graphite and diamond on substrates. 
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Table 1 . Summary of some possible gas phase reactions and energy exchange mechanisms for hydrogen 
microwave plasma 


No. 

Type 

Desig- 

nation 

Reaction Equation 

Used in 
terms 

Ref. 

1 

electron- 

translation 

mm 

e + H 2 — » e + H 2 

Qe-T, ©L 

12 

2 

electron- 

vibration 

e-V 

e + H 2 (v) -» e + H 2 (w) 

Qe-v> 0L 

12 

3 

electron- 
elec tronic- 
vibration 

E-V 

e + H 2 (v) -> e + H 2 (B 1 Z + u , C^u ) 
-> e + H 2 (v) +hv 

C^vQtad 

0L 

12 

■ 

electron- 

dissociation 

e-D 

e + H 2 (v) — > e + B- 2 * — » e + 2H 

k f. Qd 

kd, ©L 

12 

5 

dissoc.excit. of 
atom by electron 

e-D 

e + H 2 (v) -> e + H 2 * -> e + H + H* 

kf, ©l 
kd* 

15-17 

Janev 

2.2.6 

6 

ionization of 
molecule by 
electron 

e-Im 

e + H 2 (v) — > c + H 2 + + c 

kim* ©L 

■ 

■ 


e-Ia 

e + H — > e + H + + c 

kia> 

15 

■ 

excitation of 
atom 

excitation 

e + H -» e + H*(n=3) 

©L Qrad 

15 

9 

dissociative 

recombination 

drv 

e + H 2 + (v) H + H* 

kdrv 

23 

10 

vibration- 

translation 

coupling 

■ 

H 2 (v) + H 2 -> H 2 (v-1) + H 2 

Qr-v 

11 

11 

vib. quenching 
by atoms 


H 2 (v) + H -> H 2 (w) + H 

Qr-v 

2 

12 

3 body recomb, 
of atoms 

recomb 

2H + M — > H 2 + M M=H 2 ,H 

kb 

24 

13 

3 body recomb, 
of atom ions 

deioniza- 

tion 

2e + H+ -> e + H 

krl 

18, 20 

14 

radiative recomb, 
atom 

deioniza- 

tion 

e + H+ hv + H 

knrl 

18,20 

15 

dissociative 

attachment 

da 

e + H 2 -> H* + H* 

kda 

27 

16 

charge exchange 

cx 

H + H 2 + — > H 2 + H+ 


22 

17 

exc. prod.H 3 + 

e23 

H 2 + + H 2 -> H 3 + + H 

ke23 

22 

18 

vib. exc. by H + 


H+ + H 2 (v=0) -4 H+ + H 2 (v>0) 

3.2.2 

15 


































































































Table 2 - Formulas for relaxation times and collision frequencies for a moderate pressure hydrogen plasma* 
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tExcept as noted calculated from cross sections in Buckman and Phelps 12 
§Gorse, et al. 3 









Table 3 - Reaction rates for a moderate pressure hydrogen plasma 


Reaction 

Constants 
k = ATne' T a^ r 

Temp. 

Range 

Gov. 

Temp 

Source of 
Data ! 


A, 

cm 3 /s/mole 

T\ 

T a .K 




e + H 2 (v) e + H2 + + e 
kirn 

1.18xl0 16 

0 

191500 

<10 4 - 5X10 4 K 

T e 

from cross 
sections 15 

e + H-»e + H + + e 

kia 

l.OxlO 16 

0 

157583 

clO 4 - 5X10 4 K 

T e 

from cross 
sections 15 

e + H 2 (v) e + H 2 * e + 2H 
kd 

1.26xl0 16 

0 

135694 

clO 4 - 5X10 4 K 

T e 

from cross 
sections 12 

e + H3+ 3H 

-> H2*(v> 5) + H*(n=2) 

8.0xl0 17 

3.2xl0 17 

2.5:1 

-0.404 

0 

0.1-1 eV 

T e 

Janev 15 from 
cross sections 
2.2.15a, b 

e + H 3 + -» e + H + + 2H 

1.22xl0 17 

0 

179380 

5000-50000 K 

T e 

Janev 15 from 
curve fit 2.2.16 

e + H 2 (v) -> H + H~ k da 
v> 4 

2.24X10 22 

-1.45 

9592 

5000-50000 K 

Te 

Janev 15 2.2.171 

e + H 2 + -> e + H + + H(n=2) 

1.43xl0 16 

0 

195820 

5000-50000 K 

T e 

Janev 15 2.2.13 

e + H 2 + -> H + H(n) 

9.44xl0 18 

-0.604 

0 

5000-50000 K 

T e 

Janev 15 2.2.14 

e + H2 + -» e + H + + H 

1.46xl0 17 

0 

37460 

5000-50000 K 

T e 

Janev 15 2.2.12 

e + H 2 + (v) — ► H(ls) + H(nl) 
kdrv 

v=0 8.78xl0- 8 
v=l 1.81xl0' 6 
v=2 2.80xl0- 7 
v=3 8.64xl0’ 8 
v=®4 1.13x10"^ 
cm3/s/atom 

-0.29 

-0.50 

-0.66 

-0.32 

-0.77 

0 

0 

0 

0 

0 

0.02- 1.0 eV 

Te 

Nakashima, et 
al. 23 

2e + H + -»e + H k r i 

3.63X10 37 

cmfys/mole 2 

-4.0 

0 

300-5500 K 

T e 

Johnson and 
Hinnov 2 ^* 21 

e + H' -» 2e + H 

1.34xl0 13 

0.9 

22700. 

5000-50000 K 


Janev 15 7.1.1 

e + H + — » hv + H krrl 

1.46xl0 14 

-0.699 

0 

250-64000 K 

T e 

Massey 

2H + H 2 -> H 2 + H 2 kb2 

l.OxlO 17 

cm^/s/mole 2 

-0.6 

0 

50-5000 K 

T 

Cohen & 
Westbutg 24 

2H + H-»H2 + H kbl 

3.2xl0 15 

cm^/s/mole 2 

0 

0 

50-5000 K 

T 

Cohen & 
Westburg 24 

H2 + + H->H+ + H2 kcx 

3.85xl0 14 

0 

6 

not given 

T 

Karpas et al. 22 

H2 + + H 2 ->H3 + + H 

1.27xl0 15 

0 

0 

not given 

T 

KarpasetaL 22 

H 2 + H 2 H + H + H 2 

8.61xl0 17 

-0.7 

52530 

600-5000 K 

Tt 

Cohen & 
Westburg 24 

H2 + H-»H + H + H 

2.7xl0 16 

-0.1 

52530 

600-5000 K 

T+. 

Cohen & 
Westburg 24 

H + H“ -» e + 2H 

4.5x10* 

1.5 

698 

1 100-11000 K 

T 

Janev 13 7.3.2a 

H+H- -*e + H 2 (v) 

1.43xl0 15 

-0.146 

815. 

1100-11000 K 

T 

Janev 1 ^ 7.3.2b 

H + + H 2 -» H + H 2 + 

7.63xl0 13 

0 

32222. 

1100-11000 K 

T 

Janev 15 3.2.3 

H + + H' — > H(n=3) + Hfls) 

1.78 x 10 17 

0 

1768. 

1100-11000 K 

T 

Janev 15 7.2.3 


t These reactions actually depend on T v as well as T, but I have not found cross section data that depend on vibrational 
excitation state that are needed to calculate the vibrational temperature dependence. 

1 Rate must be reduced by Boltzman factor for 4. 
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electron distribution function. 
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Fig. 6 Hydrogen molecule ionization rate coefficient 
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Fig. 7 Recombination rates for H + + e reactions 
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